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Abstract. The exact one-to-one mapping between (spinless) Jordan- Wigner 
lattice fermions and (spin-1/2) spinous is established for all eigenstates of the 
one-dimensioual s = 1/2 XX model ou a lattice with an even or odd number 
N of lattice sites and periodic boundary conditions. Exact product formulas for 
the transition rates derived via Bethe ansatz are used to calculate asymptotic 
expressions of the 2-spinon and 4-spinon parts (for large even N) as well as of 
the 1-spinon and 3-spinon parts (for large odd A^) of the dynamic spin structure 
factors. The observability of these spectral contributions is assessed for finite N 
and for N ■—> oo. 



1. Introduction 

The exact solution of the XX model for exchange-coupled electron spins on a one- 
dimensional lattice was first demonstrated more than 40 years ago j2 IS- Among 
all integrable quantum many-body systems, the XX model is perhaps the one 
investigated most thoroughly [ni|l|SlEl|71|HlEl[ini[ni[I21[Il[ni. 

The XX interaction between nearest-neighbor spins with s = 1/2 along a chain 
of N sites is described by the Hamiltonian 

N 

Hxx = Y.'^StSf^, + Spy_,,]. (1.1) 

1=1 

The exact Jordan- Wigner mapping of this system onto a system of free, spinless lattice 
fermions has been at the root of most advances reported for this model. The many 
nontrivial properties of the XX model are accounted for by the non-local functional 
relation between spin operators and fermion operators. 

In spite of steady progress on the XX model reported since 1961, some questions 
of interest have resisted satisfactory answers. The objective of this paper is to 
communicate significant advances on two such issues: (i) We establish an exact one- 
to-one mapping between the fermion composition and the spinon composition of all 
XX eigenstates, thus linking the computationally convenient fermion quasiparticles 
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to the physically relevant spinon quasiparticles on a rigorous basis, (ii) We use the 
exact product formulas for transition rates previously derived in the framework of 
the Bethe ansatz to calculate exact and asymptotic expressions of the m-spinon 
dynamic spin structure factors at T = with m = 2,4 for even TV and m = 1,3 
for odd TV. The nature of the results raises questions regarding the observability of 
specific sets of m-spinon excitations and suggests that we distinguish a mesoscopic 
regime (TV ~ 10^ — 10"^) from a macroscopic regime (TV — > oo). 

The exclusion statistics of the fermion and spinon quasiparticles as well as the 
mapping between them is described in Sec. |21 The properties of m-spinon dynamic 
structure factors are discussed in Sec. 13 and the observability of m-spinon excitations 
in Sec. H 



2. Fermions versus spinons 



It is desirable to carry out calculations for the XX model using the fermions because 
they are free and to interpret all spectral properties in terms of the (interacting) 
spinons because the spinon vacuum, unlike the fermion vacuum, is at the bottom of 
the spectrum. Hence to calculate the m-spinon parts of the T = dynamic spin 
structure factors from exact transition rate formulas based on fermion momenta, we 
must understand in detail how the fermion and spinon quasiparticles are related. 

We consider the XX model for even TV or odd TV and with periodic boundary 
conditions. The Jordan- Wigner transform followed by a Fourier transform to the 
reciprocal lattice converts Hamiltonian 11. 1|) into 

'Hxx = cos kj c\^ Ckj ■ (2.1) 

Here c^, are fermion creation and annihilation operators. The sum {ki} is over the 
allowed fermion momenta: 

fc --m- m G I ^^'^'---'^^"^^ (evenTV/) 
-^m„ I {o^2,...,2Ar-2} (odd TV/) ' 

The number of fermions varies over the range < Nf < TV. 



2.1. Exclusion statistics 

The number of eigenstates composed of Nf fermions is described by the combinatorial 
formula 

Wf{Nf)=^'^ + ^y^ y d = N-Nf + l. (2.3) 

The number of (spinless) fermions present in an eigenstate is related to the quantum 
number (z-component of the total spin) of that state as follows: 

Nf = ^- S^. (2.4) 

The ground state of Hxx is non-degenerate for even TV. It contains TV/ — TV/2 
fermions. This state is reconfigured as the vacuum for spinon quasiparticles. The 
(spin-1/2) spinons obey semionic exclusion statistics. The number of eigenstates 
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composed of n+ spinons with spin up and n_ spinons with spin down is described 
by Haldane's generahzation of H2.3|l to fractional statistics 

W^K,n_)= n ( "^^^^r' (2-5) 

d. = ^(A^ + l)-^E("-'-'5..')- (2.6) 

The total number of spinons, 

n+ + n^^N,, (2.7) 

is restricted to the values Ng = 0,2, ... ,N for even N and to the values Ng = 
1, 3, . . . , for odd N. Eigenstates with a fixed number of spinons may contain different 
numbers of fermions and vice versa. However, from H2.4|l we infer 

n+ - = N - 2Nf. (2.8) 

Knowledge of and for a given eigenstate is thus equivalent to knowledge of Ng 
and Nf. 

In Tabled we list the numbers of eigenstates for iV = 6 (left) and N — 7 (right) 
with given Ng in the spinon representation and given Nf in the fermion representation. 
Summing the entries over Ng yields the number of eigenstates with Nf fermions, 
Eq. (|2.3|l rewritten as 

WfiNf) = ( (2.9) 
Summing the same entries over Nf yields the number of eigenstates with Ng spinons, 

WgiNg)=(^ ^7^^ )' (2.10) 
for a total number 

J2Wf{Nf) = Y.^giNg)^2^. (2.11) 

N, N, 



Table 1. Number VF(n_|_,n_) of eigenstates with Ns = n+ + n_ spinons and 
(in a different representation) A^^ = N/2 — (n+ — n_)/2 fermions for Tixx with 
N = 6 (left) and N = 7 (right). Subtotals FIM and ITTni are found in the last 
column and row, respectively. 
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The fermion vacuum is the entry on the upper right corner in both arrays. In the 
array for iV = 6 the entry on the far left is the spinon vacuum, the non-degenerate 
ground state for even N. The two entries on the far left in the array for N = 7 are 
1-spinon states. The fourfold degenerate ground state for odd N consists of two states 
from each entry. 



2.2. Spinon momenta 

The set of allowed spinon momentum values was first determined in the context of 
the Haldane-Shastry (HS) model for even N. Here we adapt those findings to 
the XX model and extend them to lattices with odd N. For this purpose we only 
use symmetries common to the HS and XX models: translations along the (periodic) 
chain, rotations about the z-axis and the reflection z ^ —z in spin space. 

The momentum states available for occupation by spinous are equally spaced at 
Ak = 2tt/N and their number is iVo = (TV — Ns)/2 + 1. Hence the (generalized) Pauli 
principle for spinous is semionic ATVq = —gANg with g = ^. For even N the 
allowed spinon momentum values are 

«» = ^"i„ = ^,^ + 2,...,7V-^. (2.12) 

Each available momentum state may be occupied by spinous of either spin orientation 
without further restrictions. The quantum numbers k and of all eigenstates are 
obtained from the spinon momenta Hj and the spinon spins aj via 

k = ko + ^Kj, S^^^ffj, (2.13) 

3 j 

where fco = for even N/2 and ko = tt for odd N/2. For odd N the allowed spinon 
momentum values are 

TT N N, N N, ^ 3N Ns 

-. = ^-. m^ = - + -,- + - + 2,. (2.14) 

The quantum numbers k and are again determined by H2.13|l but now with fcp = 
if {N + Ns)/2 is even and fco = tt if (A^ + Ns)/2 is odd. 



2.3. Fermion- spinon mapping 

To keep track of the (physically relevant) spinons in the (computationally convenient) 
fermion representation we consider a system of = 4 sites for illustration. The 
allowed fermion momenta for Nf = 0, 1, 2, 3, 4 are shown in Fig. ^a) and the allowed 
spinon momenta for Ns ~ 0,2,4 in Fig. H^b). An expanded version of Fig. H^a) is 
shown in Fig. |21 with all 2^ = 16 distinct fermion configurations with increasing Nj. 
The V-shaped line of Fig. Q^a) becomes the forked line in Fig. [3 

The exact spinon configuration is encoded in the fermion configuration as 
described in the following: (i) Consider the V or the fork as dividing the fermion 
momentum space into two domains, the inside and the outside. The outside domain 
wraps around at the extremes {fhi = A^modA^ = 0). (ii) Every fermionic hole 
(open circle) inside represents a spin-up spinon (square surrounding open circle) and 
every fermionic particle (full circle) outside represents a spin-down spinon (square 
surrounding full circle), (iii) Any number of adjacent spinons in the representation 
of Fig. 121 are in the same momentum state of Fig. H^b). Two spin-up (spin-down) 
spinons that are separated by € consecutive fermionic particles (holes) have momenta 
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Figure 1 . (a) Fermion momentum states available to Nf fermions and (b) spinon 
momentum states available to Ng spinous in XX eigenstates for N = 4. All 
momenta are in units of tt/N. Fermion momentum states can be either va<;aiit 
(open circle) or singly occupied (full circle). The particular fermion configuration 
shown in each row represents the lowest energy state for given Nf. Spinon 
momentum states can be either vacant or occupied by up to Ng spinous with 
arbitrary spin polarization. No specific spinon configuration is shown. 
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Figure 2. Fermion configurations of all eigenstates for N = A. Fermionic particles 
(holes) arc denoted by full (open) circles. Spinon particles with spin up (down) 
are denoted by squares around open (full) circles. The fermion momenta m; (in 
units of Tv/N) can be read off the diagram. The spinon momenta m; (also in units 
of tt/AT) and the spinon spins crj are given explicitly and can be inferred from the 
fermion configuration as explained in the text. Also given are the wave number 
k (in units of 2'k/N) and the energy E of each eigenstate. 



separated by iir /N . (iv) The spinon momenta are sorted in increasing order from the 
right-hand prong of the fork toward the left in the inside domain and toward the right 
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with wrap-around through the outside domain. 

These rules determine the exact spinon spin and momentum configuration for 
the given fermion momentum configuration of all eigenstates pertaining to an XX 
chain with A'^ = 4. Generalizations to other even values of TV are straightforward. In 
applications of these rules to odd N we must take note of the different set of allowed 
spinon momentum values and of the fact that the number of spinous is odd. 



3. Dynamic spin structure factors 

The dynamic spin structure factors for a finite system (with even or odd N) can be 
expressed as a weighted average 

S,„,{q,(^) ^ ^^e^'^^"S^,i,{q,uj)n, fJ.i' = -\ — , — \-,zz (3.1) 

n 

of functions 

S^,v{q,^)n = 27r ^ Af,^„, (q)(5 (w - w„„,) , (3.2) 

n' 

where Z ^ ^-P^n t^ii]^ p ~ (kBT)~^ is the canonical partition function and 

MLM)^\{n\S^,\n')? (3.3) 
are transition rates for spin fluctuation operators 

1 ^ 

^." = 7^E^^''^" (3-4) 
^ 1=1 

with q = (27r/iV)m, to = 0, 1, . . . , TV - 1. 



3.1. Transition rates 

Given the fermion momentum configurations 

{^1 I ■ ■ ■ J ^_r}: {^1 , ■ . . ,kj. } (3-5) 

of the initial state {n\ and final state \n') involved in 1)3. 3|l . the transition rates have 
the following explicit form 

(fisin^Mfflsin^^S^) 

ML'il) = „ : , , ^-Sr,R+lS,.Q, (3.6) 



where 



i=i j=i 

M-^,{q) = ^ ; ^-SrM^lS,,Q, (3.7) 



R 

n 

i=i »=i 



UN^U sin' 



ML'{q) = ^Sr,RS,,Q, (3.8) 



r R 

J2 ~ ^1 ^TT (3.9) 

i=\ j=i 
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is the momentum transfer during the transition, and where u„„' = {S^Y with 
— N/2 — R from H2.4|l if the two sets of fermion momenta (|3.5|l are identical, 
Unn' = 1 if they differ by exactly one element, and u„„' = if they differ by more 
than one element. The energy transfer is 



The result H3.8|l is elementary in the fermion representation. The product formulas 
H3.t)|l and (|3.7|l were derived in Ref. |T2| from determinantal expressions obtained via 
algebraic Bethe ansatz for the planar XXZ--moAe\ In the Bethe ansatz 

context the A:" and fc" were magnon momenta, i.e. solutions of the Bethe ansatz 
equations. The general solution was shown to produce (i) real, non-critical magnon 
momenta, identical to the fermion momenta p. 211 . (ii) real or complex critical pairs 
of magnon momenta (with fcj + kj = tt), different from any fermion momentum in the 
set H2.2|l . Critical pairs only occur for even N . 

The determinantal expressions from which the product formulas (|3.5|l and H3.7|l 
were derived contain factors cosi(fci + kj) in the denominator, which are zero for 
critical pairs. For this reason the conversion into the product expressions as reported 
in Ref. proceeded with the restriction that no critical pairs were present. Since 
the product expressions do no longer contain any such vanishing factors, it is very well 
possible that their validity is broader in scope than the original derivation suggested. 
On the basis of numerical tests including the sum rule 



we are indeed led to conjecture that the product expressions (|3.6|l and H3.7|) are 
universally valid for even or odd N provided we replace the critical magnon momenta 
by the corresponding fermion momenta. A proof of this conjecture will either have to 
be an extension of the existing derivation within the framework of the Bethe ansatz to 
include the critical magnon momenta or it will have to be an independent derivation 
within the fermion representation. 

3. 2. Transitions from ground state for even N 

AtT — Q only transitions from the lowest energy level survive in expression H3.1|l . For 
even N the lowest level is the (non-degenerate) spinon vacuum. Its fermion momentum 
configuration for N = 12, now denoted by the set . . . , ^^^2}' illustrated in row 
1 of Fig. 13 From that state the operator and S~ reach different sets of excitations. 
Without loss of generality we assume that N /A is an integer. 

In Szz{<l,^)o the only non- vanishing transition rates H3.8|l pertain to the set of 
2-spinon excitations (with ri+ = n_ = 1). In the fermion representation they are 
1-particle-l-hole states relative to the ground state, generated from row 1 of Fig.|21by 
moving exactly one fermion from the inside domain to the outside domain. The 
total number of such states is iV^/4, each contributing the amount 1/A^ to the 
spectral weight. Parenthetically, let us mention the polymer fluctuation operator 
P^*^,M = 1,2,3,... introduced in Ref. j^U). This operator only allows transitions 
from the spinon vacuum to the set of 2M-spinon excitations (with n+ = n_ = Af), 
i.e. to the set of Af-particle-M-hole excitations in the fermion representation. 



r 



R 




(3.10) 




(3.11) 
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Figure 3. System with N = 12 sites. Fermion momenta ki = (iT/N)rhi for tlie 
lowest energy states at Nf = N/2 (row 1) and Nf = N/2 — 1 (row 2). Also 
shown are the fermion momenta for a generic 2-spinon state (row 3) and a generic 
4-spinon state in the subspace with Nf = N/2 — 1. Also shown are segments of 
the two-pronged fork described in the context of Fig. |21 

Our main focus here is on S* ^(g,w)o. Note that we have 5 ^-{q,uJ)o — 5_| {q,(^)o 

due to reflection symmetry in spin space, even though the two functions involve 
different sets of excited states. Afl excitations in the invariant subspace with Nf = 

N/2 — 1 contribute spectral weight to S ^{q,uj)o. For their systematic generation it 

is useful to start from the state with the lowest energy in that subspace. It has spinon 
composition n+ — 2,n^ — and its fermion momenta {fcj } are illustrated in row 2 
of Fig. El 

The complete set of 2-spinon states in the same invariant subspace can now be 
described by two even- valued integer parameters N/2 < /ii < /i2 < 3N/2 as illustrated 
in row 3 of Fig. |21 They mark the two vacancies in the array of fermion momenta 
{k^} in the inside domain. These vacancies represent spinous with spin up (ri+ = 2). 
The spinon momenta are determined by the rules given in Sec. [3 

The 4-spinon excitations relevant for S ^.{q,uJ)Q have — 3 and n_ = 1. 

The spin-up spinous are described by three even- valued integer parameters N/2 < 
(Ui < ^2 < < 3N/2 as shown in row 4 of Fig. El The spin-down spinon is 
described by the even-valued integer parameter /x* in the outside domain i.e. with 
range < /i^ < N/2 or 3A^/2 < fil < 2N. In generalization to this recipe, the 

m-spinon excitations contributing to S \^{q,Lo)o have n+ = m/2 -\- = m/2 — 1, 

generated by moving m/2 — 1 fermions from the inside domain to the outside domain 
(with parameters /i*, . . . , /i^^2--i) ^^'^ thus leaving m/2 + 1 vacancies inside (with 
parameters ^i, . . . , ^J.m/2+l)■ 

3.3. Total m-spinon intensity for even N 

How is the total spectral intensity in the dynamic spin structure factor S \-{q,Lo) 

divided among the m-spinon excitations? From a previous study based on algebraic 
analysis PI] we know that in the axial regime of the XX Z antiferromagnet most of 
the spectral weight is carried by the 2-spinon excitations alone. On the other hand, 
previous Bethe ansatz studies |EI yielded evidence that in the planar regime the 
2-spinon contribution approaches zero in the limit iV — > oo. 

By using the transition rate expressions (|3.7|l we have produced numerical data 
for the relative 2-spinon intensity 12/Itot and the relative 4-spinon intensity I^/Itoty 
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Figure 4. Relative overall intensity of the 2-spinon excitations (o) and 4-spinon 
excitations (•) in S |-(5,cj)o plotted versus l/N. 



where 

q n m— 2,4,... 

The data are shown in Fig. 21 We see that the 2-spinon intensity accounts for more 
than half the total intensity in chains with up to at least N ~ 300 sites. Nevertheless, 
the data are consistent with the conclusion that the relative 2-spinon intensity drops 
to zero in the limit TV — > oo. 

The relative contribution of the 4-spinon excitations first rises with N. It levels 
off at ~ 48% for N ~ 2000 and then starts to decrease. In all likelihood, the 
relative 4-spinon intensity and, for that matter, any relative m-spinon intensity taken 
individually, will vanish in the limit N oo. However, the combined 2-spinon and 
4-spinon relative intensity stays above 80% for chains with up to ~ 2000 sites. 



3.4- Asymptotic transition rates for even N 

How is the m-spinon intensity of S y {q, lu)q distributed in the (g, CL;)-plane for large iV? 

To answer this question we must determine the explicit dependence of the transition 
rate expression (|3.7I) on the spectral parameters q and uj. The product nature of 
H3.7(l suggests the following factorization of the transition rate expression between the 
spinon vacuum and an m-spinon state (for m = 2, 4, . . .): 

^ « ' ^ . ' 

scaled transition rate scaling factor 

Here {fc^*} and {k^} are the fermion momenta of the lowest state in the invariant 
Hilbert subspaces with N/2 fermions and N/2 — 1 fermions, respectively. The fermion 
momenta {fc™} of an arbitrary m-spinon excitation differ from the set {k^} as 
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explained in the context of Fig. |31 In the representation H3.13|l we use the transition 
rate between the two reference states (second factor) as a common scaling factor for all 
m-spinon excitations. In the scaled transition rates (first factor), huge cancellations 
take place because different m-spinon excitations differ by no more than m fermion 
momenta. 

The asymptotic scaling factor must be calculated by expansions of the product 
expression 

M-({fcO}, [kr]) = ^Cn{N/2), (3.14) 

sin^'-S(^ ^ 1/2) . TT 

CN{n)=YY . 4;_i ■ , "n^iTj (3.15) 

sm -ql N 

as inferred from H3.7|l for the two specific states. The result is of the form 

M-({fc«},{fc,r}) ^Va^C, (3.16) 

where here and henceforth the symbol denotes an asymptotic equality for large 
N that ignores corrections of the kind x[l + 0{N^^)]. The constant 

C= lim Cn{N/2) (3.17) 



oo 



= V^exp + 6C'(-1)^ 0.73739071 . . . , 

where C,{z) is the Riemann zeta function is familiar from previous calculations of 
correlation functions for the XX model by a different approach We now present 
explicit results for the 2-spinon and 4-spinon transition rates. 

3.5. Asymptotic 2-spinon transition rates 

The scaled transition rate for a 2-spinon excitation can be reduced to the following 
product: 



M-({fcO},{fcf}) 



with the function 



= sin^ "^^-^ n ^{^mN/'i - ^o) (3.18) 



1-2/7 1\ 

sm rji 

depending on two parameters /^i and /i2 via 

u N 

V'.-y-j, =W» z = l,2. (3.20) 

Performing the limit ^ oo is delicate because of singularities along the spectral 
boundaries of the 2-spinon excitations. To leading order in N , expression (|3.18|l 
becomes 



M-({fcO},{fc-}) iV2 COSK1COSK2 
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where the two parameters Ki , K2 are related to the energy transfer to and momentum 
transfer q as follows: 

2 2 
(jj = — cos Ki, q = TT — Ki mod 2tt. (3.22) 

i=l i=l 

Note that the limits {kf} — > {k^} and N ^ 00 are not interchangeable. Hence 
we cannot recover the unit value of the scaled transition rate by taking the limit 
Ki — K2 — > TT and Ki + K2 — 27r in the asymptotic expression H3.21|l . 

Assembling the results (|XTB|) and with lIX^ into ifTT^ yields the 2 -spinon 

transition rate function 

Af-t 4C 4sin^(g/2)-^^ 

M2 (g,w) = — ^ — (3.23) 

over the (asymptotic) range | sinq| < lo < 2\ sin(g/2)| of the 2-spinon spectrum. The 

2-spinon part of the dynamic spin structure factor S ^(g,a')o is then obtained by 

multiplying the 2-spinon transition rate function l|3.23|) with the 2-spinon density of 
states, 

D2[q,u:) ^ , ^ (3.24) 

y'4sin2(g/2) -a;2 

yielding 

2C v/4sin^(g/2)-^2 
^ c^-sin^, ■ 

The interpretation of the asymptotic result (|3.25() requires caution. It is indicative of 
the 2-spinon spectral-weight distribution in the following sense: in a histogram plot of 
S_!^{q,uj)o based on the transition rates 13.7|l . the shape of (|3.25|) comes into better 
and better focus as N grows larger, while the relative intensity of this structure fades 
away. 

3.6. Asymptotic A-spinon transition rates 

We start again from the factorized expression H3.13|l . The asymptotic scaling factor 
remains the same: Eq. (|3.16|) . The surviving factors in the scaled transition rate after 
the massive cancellations now depend on four parameters (see Fig. O : 

3 

TTsin^^i^ 

M-({fc,0},{fc,f}) ..2.2(^1 A^<^ 
— N sm ' ' 



M-{{k^},{kl}) \2 4;^ 



sm 

1=1 

3 



X n - ^,) (3.26) 



where 



cbi3N/4~^il/2) f \ 



^j = y-^, K^=l^fJ,^, i = 1,2,3, K\=rinl. (3.27) 
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From this result we have derived the foUowing asymptotic expression for the scaled 
4-spinon transition rates: 

Asin^^l^ 
M-iimAkf}) a 4 cosK) 2 

_[J_cos(k,) _[J_sm 

4=1 1 = 1 

The asymptotic dependence of the energy and momentum transfer of 4-spinon 
excitations on the four parameters is 

3 3 

uj = COS kI — cos Ki, cj ~ TT + kI — Ki mod 27r. (3.29) 

1=1 i=l 

From inspection of the 2-spinon result (|3.21|) and the 4-spinon result H3.28|l we can 
already see features of the general structure for the m-spinon result emerging. 




Figure 5. Excitation energy versus wave number of all 2-spinon (o) and 4-spinon 
(■) excitations for A'^ = 64. Tiie solid lines are the spectral boundaries for A'^ — > oo. 



The range of the 4-spinon excitations in the {q, w)-plane in relation to the range 
of the 2-spinon excitations is illustrated in Fig. |S1 for TV = 64 and asymptotically for 
— > oo. We can discern the nonuniform density of 2-spinon states H3.24|l in the 
distribution of circles. Furthermore, from the distribution of dots we see that the 
density of 4-spinon excitations is very small near the spectral threshold and remains 
relatively small at energies uj < 1. In the limit N ^ oo, the 2-spinon and 4-spinon 
excitations have a common spectral threshold, 

e2L(g) = e4L(g) = I sin^l, (3.30) 
and different upper boundaries. 



. q 

sm — 
2 



e2uiq) = 2 
£4(7 ((7) = 4 max 





. 9 




. q - 27r 






sm — 




sm 






4 




4 





(3.31) 
(3.32) 
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3. 7. Transitions from the ground state for odd N 

For odd N the ground state is fourfold degenerate. Two of its vectors (labeled A and 
A') are located in the invariant subspace with Nf = {N — l)/2 fermions and the other 
two vectors (labeled B and B') in the subspace with Nf = {N + l)/2 fermions. We 
then have to consider transitions from each of the four ground-state vectors to the 
sets of 771-spinon excitations with m = 1,3, ... in the accessible invariant subspaces. 
Without loss of generality we assume that {N + l)/4 is an integer. 

The fermion momentum configurations of two of the four ground states are 
shown in rows 1 and 4 of Fig. for iV = 11. It will suffice to consider these 
two vectors (with fermion momenta {/cj, . . . , ^^'^ i^i' ■ ■ ■ ' ^^^d 

calculate from them the functions Sfj,iy{q,uj)A and Sp,v{q,uj)B- Symmetry dictates 
that the other two ground state vectors yield the fimctions Sfj,u{q, (-o)a' — Sfj,u{^q, lli)a 
and S^ij{q,ijj)B' = Sfiu(—q,^)B- Hence the T — Q dynamic spin structure factor l|3.1|l 
for odd N can be rewritten as 

Sf_,^{q,Lu) = i [Sf^i,{q,uj)A + Sf,^{-q,uj)A + S^^{q,uj)B + S^„{-q,Lo)B] ■ (3.33) 
Different sets of excitations are again relevant for pLv — zz,~\ — , — h. 
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Figure 6. XX chain with A'^ = 11 sites. Fermion momenta ki = {■K/N)mi for 
one of two lowest energy states at Nf = (iV — 1)/2 (rows 1 and 5), Nf = {N — 3)/2 
(row 2), and Nf = {N + l)/2 (row 4). The other state in each case (henceforth 
referred to by subscripts A' or B') is obtained via reflection at the line rhi = 11. 
Also shown are the fermion momenta for a generic 1-spinon state in the subspace 
with Nf = {N — l)/2 (row 6), a generic 3-spinon state in the subspace with 
Nf = {N — 3)/2 (row 3), and a generic 3-spinon state in the subspace with 
Nf = {N - l)/2 (row 7). 

In Szz{q, ^^)a the non- vanishing transition rates H3.8II now include all {N+l)/2 1- 
spinon states (with n_|_ = 1), including the ground-state vector A itself, and a subset of 
{N — 1)^/4 3-spinon states (with — 2,n^ — 1), all in the same invariant subspace. 
In row 1 of Fig. |H| the 1-spinon states are obtained by moving the vacancy to any 
site in the inside domain, whereas the contributing 3-spinon excitations are obtained 
by leaving the vacancy inside as is and moving one particle from anywhere inside to 
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0.02 



Figure 7. Relative overall intensity of the 1-spinon excitations (I) and 3-spinon 
excitations (o) in S- + {q, u))b > a-^d of the 3-spinon excitations (•) in S--\-{q, 
all plotted versus 1/A'^. 



anywhere outside. The relevant spectrum of Szz{.q,^)B is identified in hke manner by 
starting from row 4 and interchanging the roles of the two domains and the roles of 
particles and vacancies. 

As in Sec. 13.21 our main focus is on = — h. The spectrum for S \-{q^Lii!)A 

includes all eigenstates with Nf = {N — 3)/2 and the spectrum for S ^{q,uj)B all 

eigenstates with A^/ = (A^— 1)/2. For a systematic generation of all relevant excitations 
we again consider one state each in the invariant subspaces with Nf = {N — 3)/2 and 
Nf ~ {N — l)/2 as reference states for the scaled matrix elements. These two states 
have fermion momentum configurations {k^ , . . . , ^(7v-3)/2}^ ^^"^ i^i ^ • ■ • > ^(^n-i)/'^^b 
as illustrated (for A^ = 11) in rows 2 and 5, respectively, of Fig.|S| 

Reference state is a 3-spinon state (with n+ — 2, n_ — 1) whereas reference 

state is a 1-spinon state (with n_|_ = 1). The remaining states in the subspace 

of {k~}A are m-spinon states for m = 3, 5, . . . with n+ = (m + 3)/2, n- = (m — 3)/2, 
whereas the remaining states in the subspace of {k~}B are m-spinon states for 
m = 1,3,... with = {m + l)/2,n_ = (m — l)/2. The 1-spinon and 3-spinon 
states thus generated are illustrated for A^ = 11 in rows 3, 6, and 7 of Fig. El Note in 
particular that 1-spinon excitations only occur in S — |-(g, a;)^ but not in 5* — \.{q,uj)A- 
We have produced extensive finite-A^ data for the relative 1-spinon intensity 
^1 I ^tot ^'^d relative 3-spinon intensities j and /Itot i where 

N - 1 



q n 



771 — 3,5,. 

E 

rn— 1,3, . 



2 ' 
N + 1 



(3.34) 



The data are shown in Fig. We see that the contribution of the 1-spinon 
excitations to the total intensity remains significant for chains with up to several 
hundred sites and the contribution of 3-spinon excitations for chains up to several 
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thousand sites. Nevertheless, the data give a clear indication that both contributions 
vanish in the limit N oo. 

3.8. Asymptotic transition rates for odd N 

In the factorization H3.13|l we now have to distinguish two scaling factors for the 
degenerate ground state. Both must be calculated anew, from the product expressions 

Af-({fc,"}A,{fcr}^) ^ ViVCAr((iV-l)/2)cos2|, (3.35) 



Af-({fcO}B,{fcr}s) = ViVCAr((iV-l)/2), (3.36) 
with Cn{(N — l)/2) defined in (|3.15() . The asymptotic results are again of the form 

M-{{k^}AAK}A) - M-{{knB,{K}B) = VnC (3.37) 

with the constant C from H3.17II . 

The scaled transition rate for a 1-spinon excitation is reduced to 

7Vcos2(ryi/;i)0(V'i)0((7V - l)/2 - Vi), (3.38) 



M-{{mB,{k}}B) 2, 



M-{{mB,{K}B 

where 



ipi = — — , Ki = rifii. (3.39) 



2 4 

Asymptotically for N ^ oo Eq. (|3.38|) becomes 

M-{{k^}BAkl}B) a 1 ^ 
= -TT cot 



M-({k°}B,{K}B) N 

The 1-spinon energy and momentum transfers are 

a a 

(jj = — cosKi, ~ ~2 — (3-41) 

Note again that the limits {kl}B {^^Jb and N ^ oo are no longer interchangeable 
in the asymptotic result. The 1-spinon transition rate function thus becomes 

Mr(g,t^)B = ^tan|. (3.42) 

Taking into account the fourfold ground-state degeneracy, we can write the asymptotic 
expression for the 1-spinon part of the dynamic spin structure factor in the form 

S^:}l{q,uj) = ^^tan|5(cj- Isinql). (3.43) 

Different sets of 3-spinon states are reached via from the ground-state 
components A and B. Hence we must work with two different results for the scaled 
transition rates. For ground-state B we obtain the expression 

M-mBAkf}B) =^tan^(,^*),i,2(^(^^_^^))-Q_ cos^W,) 



'^(-V'l*) 



cbiiN~l)/2-ri) 



^Y[mmN-l)/2~yj,) (3.44) 
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depending on the three parameters 
/ii iV + 1 



(3.45) 



2 4 

The asymptotic expansion of (|3.44l) yields the foUowing leading term, now 
parameterized by ni — rjfii, K2 — ??M2: a-nd — ri^\: 



. 1 sm (^^) tan(---) 



(3.46) 



2 4^ 

i=l i=l 



The energy transfer and momentum transfer in the asymptotic regime are 



= cos Kl — COS Ki, q = k\ — Ki mod 27r. (3.47) 



i=l 




Figure 8. Excitation energy versus wave number of all 1-spinon excitations (o) 
and 3-spinon excitations (•) from ground-state B for = 63. The solid lines are 
the 3-spinon continuum boundaries for Af — > oo. The corresponding 1-spinon and 
3-spinon spectrum from ground-state B' is the mirror image reflected at q/iv = 1 
of the spectrum shown. 



The range in the (g, a;)-plane of the 1-spinon and 3-spinon excitations from 
ground-state B is shown in Fig.|Slfor N = 63. The corresponding spectrum reached 
from ground-state B' is the one with wave numbers q replaced by 2tt — q. In the limit 
N —^ oo the single branch of 1-spinon excitations coincides with the lower boundary 
of the 3-spinon continuum, 

ei(g)-e3L(g) = |sing|. (3.48) 

The upper boundary of the 3-spinon continuum is 

(3.49) 



€3(7(9) = 3 max 





. 9 




. q - 27r 




sm — 




sm 
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The scaled transition for 3-spinon transitions from ground-state A reads 

3 



i<j 



M-({fcfU,{fe°}^) _ ^ 

M-{{KU,{knA) 2 3 



7 = 1 



n<^(V',-M^+i)/2-v,) 
1=1 



depending on the three parameters 



7/>, = y z= 1,2,3. 



(3.50) 



(3.51) 



z 4 

The leading term of the asymptotic expansion becomes a function of the three variables 



ncos(---Jsm 



with associated energy and momentum transfer 



a \ ^ o 37r \ ^ 

= — ^^cosKj, ^ ~ ~2 A^'^^ ™° 

i=l i=l 



27r. 



(3.52) 



(3.53) 




Figure 9. Excitation energy versus wave number of all 3-spinon excitations (•) 
from ground-state A for N = 63. The solid lines are the 3-spinon continuum 
boundaries for A'' — » oo. The corresponding 3-spinon spectrum from the ground- 
state component A' is the mirror image reflected at g/jr = 1 of of the spectrum 
shown. 
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The range in the {q, ti;)-plane of the 3-spinon excitations from ground-state A is 
shown in Fig. O for N — 63. The corresponding spectrum reached from ground-state 
component A' is the one with wave numbers q replaced by 2TT — q. In the hmit N oo 
the 3-spinon continuum boundaries for states A and A' combined are the same as those 
of states B and B' combined. 

4. Conclusion and outlook 

The two main results on the XX model reported in this paper are (i) the exact 
one-to-one mapping between fermion and spinon compositions of all eigenstates and 
(ii) the applications of exact product expressions for the transition rates as needed 
for the calculation of dynamic spin structure factors. The mapping has provided us 
with a framework that allows us to explore any property of interest by using the 
computationally convenient (free) fermion quasiparticles and interpret the results in 
terms of the physically relevant (interacting) spinon quasiparticles. The very product 
nature of the transition rate expressions for the XX model has put us at a huge 
advantage in the face of the equivalent determinantal expressions known for the XX Z 
model in general, when it comes to the evaluation of explicit numerical results or the 
calculation of asymptotic analytic results for large systems. 

The applications worked out in Sec. |21 call for a new assessment of some issues 
and set the stage for the pursuit of new questions of interest. Chief among the issues 
that need to be reevaluated in the light of the results reported here concerns the role 
in general and the observability in particular of specific sets of spinon quasiparticles 
in dynamic spin structure factors and related quantities. 

The evidence presented in Sec. |31 suggests that the relative intensity of any 

individual m-spinon contribution to the dynamic spin structure factor S \-{q,iLj) at 

T = vanishes in the limit N ^ oo. This conclusion, which is expected to hold 
true not only for the XX model but within the entire planar regime of the XX Z 
model is in contrast to what is known on rigorous grounds for the XX Z model in 
the axial regime including the isotropic limit {XXX model): the relative intensity of 
the 2-spinon excitations alone amounts to 73% in the Heisenberg limit and gradually 
increases to 100% in the Ising limit |22II21| . 

The dominant singularities in the dynamic structure factor have been shown to be 
extractable from the 2-spinon part alone in the XXX model 1221. This is manifestly 
not the case in the XX model and unlikely the case for the planar XXZ model 
in general. A fair amount of exact information about the singularity structure of 

S ^{q, w) at T = can be gleaned from existing results important 

task will be to assemble the asymptotic m-spinon results from the approach taken here 
to recover known exact results for iV ^ oo. A hopeful sign for this endeavor is that the 
constant C as given in H3.17|l . which plays a crucial role in exact results for ^ oo, 
already emerges from the asymptotic m-spinon expressions worked out here. 

Dynamic spin structure factors are accessible more or less directly to experimental 
investigations via neutron scattering, NMR, ESR, light scattering, and other probes 
in the context of quasi-one-dimensional magnetic compounds. Even small amounts 
of impurities or other defects will turn a crystalline sample into an ensembles of 
mesoscopic chains of various lengths (A^ ^ 10^ — 10"^). Under what circumstances 
does the distinction between a physical ensemble of mesoscopic chains and a statistical 
ensemble of macroscopic chains matter? 

If the relative intensity of the 2-spinon contribution is dominant for mesoscopic 
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chains as well as for macroscopic chains, which is the case for the XXX model and the 
axial XX Z model, then the distinction has minimal impact. If, on the other hand, the 
relative 2-spinon intensity is dominant for chains of length N ~ 10^ and the relative 
4-spinon intensity is dominant for chains of length N ~ 10'^ but both contributions 
taken individually or combined become negligible in macroscopic chains, as is the case 
in the XX model, then the distinction cannot be ignored. Hence the significance of 
the results presented in Sec. 01 

In conclusion we mention two further unresolved issues to which the work reported 
here has led. They are the subject of work currently in progress. A first set of questions 
calls for a derivation of the product expressions H3.6I) - (|3.7|) and for the establishment 
of connections to the many exact results for dynamic correlation functions already 
known for the XX model. A second set of questions calls for a detailed study of the 
spinon interaction and for the exploration of the physical relevant quasiparticles whose 
physical vacuum is the XX ground state in a magnetic field. 
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